PROCEEDINGS 

OF SCIENCE 



Spectral properties of a sample of type 1 AGNs: 
influence of star formation 



Luka C. Popovic, Jelena Kovacevic and Milan S. Dimitrijevic 

Astronomical Observatory 
E-mail: 



lpopovicgaob . rs 



To find the spectral properties of AGNs in optical spectral band (around the H/3 line) we con- 
structed an Fe II template, that is covering Fe II emission from 4100 A to 5600 A. Using the new 
Fe II template we explored spectral properties of 302 type 1 AGNs. The most interesting results 
we found is the correlation of the Hj3 Full Width at Half Maximum (FWHM) and luminosity 
for a subsample of type 1 AGNs where the ratio of narrow lines indicates a significant starburst 
contribution. 



Workshop on: Nuclei of Seyfert galaxies and QSOs - Central engine & conditions of star formation, 
November 6-8, 2012 

Max-Planck-Insitut fAijr Radioastronomie (MPIfR), Bonn, Germany 



* Speaker. 

'The Fe II template can be found at: http://servo.aob.rs/jelena/link8.html 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Spectral properties of a sample of type 1 AG/VLuka C. Popovic, Jelena Kovacevic and Milan S. Dimitrijevic 



1. Introduction 

There are some indications that evolution of AGNs is probably related with starburst regions 
(see Lfpari and Terlevich 2006; Mao et al. 2009; Sani et al. 2010). Namely, it is possible that 
AGNs in an earlier phase of their evolution are composed of starburst (star-forming) regions and 
the central engine, that is AGN, and, during evolution, the starburst contribution becomes weaker 
and/or negligible (Wang and Wei 2006, 2008; Mao et al. 2009). On the other hand, it is expected 
that Fe II strength may vary with evolution of AGNs. The optical Fe II (XX 4000-5400 A) emission 
arises from numerous transitions of the complex Fe II ion. Numerous lines are overlap, forming 
the complex shapes in XX 4000-5400 A range and making one of the most interesting features in 
AGN spectra. It is seen in almost all type-1 AGN and it is especially strong in narrow-line Seyfert 
Is (NLS Is). Origin of the optical Fe II lines, the mechanisms of their excitation and location of the 
Fe II emission region in AGN, are still open questions (see Kovacevic et al. 2010 and references 
therein). There are also many correlations between Fe II emission and other AGN properties which 
need a physical explanation. 

Lfpari and Terlevich 2006 have explained some properties of AGN by an "evolutive unification 
model". In this model, accretion arises from the interaction between nuclear starbursts and the 
supermassive black hole. Thus, young AGN have strong Fe II, BALs, weak radio emission, the 
NLR is compact and faint, and broad lines are relatively narrow. In contrast with this, old AGNs 
have weak Fe II, no BALs, strong radio emission with extended radio lobes, the NLR is extended 
and bright, and the broad lines have greater velocity widths. 

In this paper we analyze how the presence/absence of the starburst contribution to the AGN 
spectra influences to the correlations between different spectral properties, specially to the Fe II 
emission. To investigate the iron lines, we use our model of iron template, based on the physical 
properties of iron region. 

2. The sample and analysis 

For this investigation we use the sample of 302 AGNs Type 1 , selected with criteria described 
in Kovacevic et al. (2010). 

The sample has approximately an uniform redshift distribution in the range of z=0-0.7 and 
negligible host galaxy contribution. 

As it is described in Kovacevic et al. (2010), the emission lines within the 4400 - 5500 A 
range are fitted with multiple Gaussians, where each Gaussian represent contribution from different 
emission region. The Balmer lines are fitted with three components: a narrow, an intermediate and 
a very broad component (H/3 NLR, ILR and VBLR, respectively). We assumed that all narrow lines 
originate from the same emission region, and consequently they have the same width and shift. The 
H/3 broad component is taken as the sum of the H/3 ILR and H/3 VBLR components, and FWHM 
(Full Width Half Maximum) of H/3 line is measured only for the H/3 broad line (see Kovacevic et 
al. 2010). The iron lines within 4400 - 5500 A range are fitted with template presented in Section 
3.1. 

The possible contribution of starbursts in AGN Type 1 spectra has been reported, using the 
BPT diagnostic (BPT - Baldwin, Philips & Terlevich 1981) in some investigations: Popovic et 
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al. (2009) found that in the case of NLSyl galaxy Mrk 493, the narrow-line ratios correspond to 
starbursts rather than to an AGN origin. Similar result is found in paper of Mao et al. (2009) for 
the three AGN Type 1, which are mainly excited from starburst in narrow lines, since they belong 
to the H II part of BPT diagram. 

Since we have complete measurements of line parameters only for the narrow H/3 and [O III] 
lines in the whole sample (302 AGNs), we accepted a criteria of R=log([0 IITJ/H/3 NLR) = 0.5 as 
an indicator of the predominant starburst emission contribution to the narrow emission lines (for 
details see Popovic & Kovacevic 2011). We divided our sample into two subsamples: R < 0.5 (91 
AGNs, hereafter starburst dominant) and R > 0.5 (210 AGNs, hereafter AGN dominant). Then, 
we performed correlations between measured spectral properties for two subsamples, in order to 
check if there are some significant differences in correlation coefficients, which may be signature 
of different physical properties of emission regions. 



3. Results 

3.1 Fe II template 

For investigation of the origin of iron lines and for analysis of their correlations with other 
lines, it is necessary to apply a good template which will fit well iron lines within AA4000-5500 A 
range. Also, precise Fe II template is needed to accurately separate the [O III] and B aimer lines, 
which overlap with Fe II. But, the construction of iron template is very difficult since the iron lines 
form the features of a complex shape. 

A few empirical and theoretical templates are proposed in literature (Veron-Cetty et al. 2004, 
Dong et al. 2008, Bruhweiler & Verner 2008), but they can not explain very well all varieties of 
the Fe II emission. Namely, in the case of some AGN spectra, specially when the blue bump of Fe 
II emission (4400-4700 A) has higher strength than red one (5100-5600 A), these models can not 
enable a good fit (see Appendix B in Kovacevic et al. 2010). 

We calculated the Fe II template, using 67 Fe II emission lines, identified as the strongest 
within the AA4100-5600 A range. The 52 of them are separated in the five line groups according 
to their lower level of transition: 3d 6 ( 3 P2)4s 4 P, 3d 6 ( 3 F2)4s 4 F, 3d 5 4s 2 6 S, 3d 6 ( 3 G)4s 4 G and 
3d 6 ( 3 H)4s 2 H (in further text 4 P, 4 F, 6 S, 4 G and 2 H group of lines). A simplified scheme of 
those transitions is shown in Fig |l[ 

The lines from five line groups describe about 85% of Fe II emission in the observed range 
(AA4000-5600 A), but about 15% of the Fe II emission can not be explained with permitted lines 
which excitation energies are close to these of lines from the five line groups. The missing parts 
are around ~4450 A, ~4630 A, ~5130 A and ~5370 A. 

There are some indications that the process of fluorescence (self-fluorescences, continuum- 
fluorescences or Lya and Ly/3 pumping) may have a role in appearing of some Fe II lines (Vemer 
et al. 1999, Hartman & Johansson 2000). They could supply enough energy for exciting the Fe II 
lines with high energy of excitation, which could be one of the explanation for emission in these 
wavelengths. 



3 



Spectral properties of a sample of type 1 AG/VLuka C. Popovic, Jelena Kovacevic and Milan S. Dimitrijevic 




4200 4400 4600 4800 5000 5200 5400 5600 
Wavelength (in A) 



Figure 1: The simplified Grotrian diagram showing the strongest Fe II transitions in the XX 4100-5600 
A region (top). Lines are separated into five groups according to the lower level of transition (middle): P 
(dotted line), F (dashed line), S (solid line), G (dash - dotted line), and H (two-dashed line). Bottom: the 
lines from the five line groups (solid line) and additional lines taken from I Zw 1 (see Kovacevic et al. 2010), 
represented with dots. 

In order to complete the template for missing 15%, we selected 15 lines, which probably arise 
with some of these mechanisms, from Kurutcz database 1 . The selected lines have wavelengths on 
missing parts, strong oscillator strength and their energy of excitation goes up to ~11 eV. Relative 
intensities of these 15 lines are addopted from I Zw 1 spectrum by making the best fit together with 
Fe II lines from three line groups. The table of selected lines and their relative intensities are given 
in Kovacevic et al 2010) and Shapovalova et al. (2012). 

We have assumed that each of lines can be represented with a Gaussian, described by width 
(W), shift (d) and intensity (I). Since all Fe II lines from the template probably originate in the 
same region, with the same kinematical properties, values of d and W are the same for all Fe II 
lines in the case of one AGN, but intensities are assumed to be different. We suppose that relative 
intensities between the lines within one line group ( 4 P, 4 F, b S, 4 G or 2 H) can be obtained as 
(Popovic et al. 2009, Kovacevic et al. 2010): 

'http://kurucz.harvard.edu/linelists.html 
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Figure 2: An example of the best fit of the AA4000-5600 A, region of Ark 564: (A) the observed spectra 
(dots) and the best fit (solid line). (B) H/3, Hy, and H5 fit with the sum of three Gaussians representing 
emission from the NLR, ILR and BLR. The [O III] AA4959, 5007 A lines are fit with two Gaussians for 
each line of the doublet and He II A 4686 A is fit with one broad and one narrow Gaussian. The Fe II template 
is denoted with a dotted line, and also represented separately in panel (C). 



I l = ( y h)A.Z±. e -(E i -E2)/kT (31) 
h A.1 fa g2 

where I\ and h are intensities of the lines with the same lower level of the transition, X\ and 
%2 are transition wavelengths, g\ and g2 are statistical weights for the upper energy level of the 
corresponding transition, and f\ and fa are oscillator strengths, E\ and £2 are energies of the upper 
level of transitions, k is the Boltzman constant and T is the excitation temperature. The list of 52 
selected lines from three line groups, their oscillator strengths and calculated relative intensities 
for different excitation temperatures are given in Kovacevic et al. (2010) and Shapovalova et al. 
(2012). 

Finally, the template of Fe II is described by 9 free parameters in fit: width, shift, six param- 
eters of intensity - for 4 P, 4 F, 6 S, 4 G and 2 H line groups and for lines with relative intensities 
obtained from I Zw 1 . The 9th parameter is excitation temperature included in the calculation of 
relative intensities within line groups. 

The fit of the Fe II template is shown in Fig |2[ We apply this template on large sample of 
AGNs (302) from SDSS database, and we found that the template fits well Fe II emission. The Fe 
II template, as well as web application for fitting Fe II lines in AGN spectra with this model, are 
given at: http://servo.aob.rs/ jelena/index.html, as a part of Serbian Virtual Observatory. 
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Table 1: Spearman rank order correlations for total sample and for subsamples: (1) log([0 III]5007/H/3 
NLR)> 0.5 (AGN dominant) and (2) log([Q III]5007/Hj3 NLR)< 0.5 (starburst dominant). 
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3.2 Influence of starbursts to spectral properties 

We performed correlations between Fe II lines and other spectral properties (continuum lu- 
minosity, broad and narrow line widths, equivalent widths of lines) for two subsamples, AGN and 
starburst dominant. We found significant differences in some correlations for these two groups of 
objects. 

We found that the width of the broad H/3 is in a significant correlation with the continuum lu- 
minosity (L5100) for the starburst dominant subsample, while there are no any correlations between 
these parameters for the AGN dominant subsample (see Fig 3). 

The relation between EW Fe II and FWHM (FWZI) H/3 is also different for the starburst 
dominant and AGN dominant subsamples (see Fig. 4, Table 1). Namely, the well known anticorre- 
lation EW Fe II and FWHM H/3, which is part of EV1 correlations from paper Boroson and Green 
1992 (BG92), is noticed only for AGN dominant subsample (r,=-0.37, P=3.1E-8), while starburst 
dominant subsample shows even opposite, but statistically not significant trend (r v =0.26, P=0.01). 
Note here that Grupe (2004) also found positive correlation between EW Fe II and FWHM H/3 for 
NLSyl objects, and the negative correlation for BLSyls. It is the same case for EW [O III] vs. EW 
Fe II anticorrelation, which is also part of BG92 EV1 correlations. Namely it is stronger in AGN 
dominant subsample (r s =-0.38, P=8E-9) than in starburst dominant (r s =-0.27, P=0.01). Fe II lines 
show significant inverse Baldwin effect (see Baldwin 1977) in the AGN dominant subsample (r s = 
0.29, P=2.4E-5), while in other subsample (r s =0.26, P=0.01) there is no correlation. 

As it can be seen in Table 1 , there are weak correlations between EW Fe II and the EW of some 
narrow lines (NLR H/3 and [O III]) only for AGN subsample. Contrary, the correlation between 
EW Fe II and EW broad H/3 is present only for starburst subsample (r s = 0.49, P=9.4E-7), while for 
AGN subsample, no correlation is found (r s =0.04, P=0.57). 

In principle, the Fe II emission is sligthly stronger (in average) in starburst dominant subsample 
than in AGN dominant. Average EW Fe II for AGN dominant subsample is: 97 A, and for starburst 
dominant: 1 19 A). 

4. Conclusions 

In this paper we performed correlations between Fe II lines and other spectral properties in 
the sample of 302 type 1 AGNs. We devided the sample in two subsamples, AGN dominant and 



6 



Spectral properties of a sample of type 1 AG/VLuka C. Popovic, Jelena Kovacevic and Milan S. Dimitrijevic 




□ log ([O lll]/Hb NLR) > 0.5 
■ log ([0 lll]/Hb NLR) < 0.5 



3.0 3.2 3.4 3.6 
log (FWHM Hp) 



rs = 0.81 
P< 1E-16 



log ([OIII]/HpNLR)<0.5 



3.2 3.4 3.6 

log (FWHM Hp) 




log ([OIII]/HpNLR)>0.5 



3.2 3.4 3.6 3.8 4.0 
log (FWHM Hp) 





■ 










° • '.n*^ a □ □ 

1 mm n aa u u □ - 

r n □□ □ 

■ □ 




□ log([OIII]/HbNLR)>0.5 
■ log([OIII]/HbNLR)<0.5 


3.8 4.0 


4.2 4.4 4.6 


log (FWZI) 



■ " rs = 0.70 

P= 1.6E-14 

log([OIII]/HpNLR)<0.5 



3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 

log (FWZI) 




□ □ □ 



log ([O lll]/H0 NLR) > 0.5 

" ^6 



4.0 4.2 4.4 

log (FWZI) 



Figure 3: Correlations between continuum luminosity and FWHM Hj3 (and FWZI Hj3) for all AGNs (panels 
left), AGNs with R <0.5 (panels in middle) and R >0.5 (panels right). 
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Figure 4: The same as in Fig. 3 but for correlations between EW Fe II and FWHM Hj3 (and FWZI Hj3). 
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starburst dominant, in order to analyze how the presence/absence of the starburst contribution to 
the AGN spectra influences to the correlations between Fe II and spectral properties. 

We found significant difference between all analyzed correlations for different subsamples, 
which reflex different physics and geometry of emission region caused by presence/absence of the 
starbursts. However, there is no big difference in Fe II emission, except that Fe II in starburst 
dominant subsample tends to be stronger. The most interesting correlation is between FWHM and 
luminosity in the subsample of type 1 AGNs where the ratio of narrow lines indicates presence of 
the star-froming regions. This may be explained by two possibilities: 1) a non conventional BLR 
or 2) that AGNs for this subsample are accreating with the similar rate. More detailed discussion 
can be found in Popovic & Kovacevic (201 1). 
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